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ABSTRACT 

We consider the emission within the fireball model framework of very high energy, ~ 1 GeV to 
> 1 TeV photons, on a minute time scale, during the onset of fireball deceleration due to interaction 
with surrounding medium. Our time dependent numerical model includes exact treatment of electron 
synchrotron emission, inverse-Compton scattering, pair production, and evolution of electromagnetic 
cascades (initiated by pair production or photo-production of pions). We find that (i) The 1 GeV- 
10 GeV flux is not sensitive to model parameters and is ~ 10~ 7 ergcm~ 2 sec -1 for z — 1 bursts, 
well within the detection capabilities of GLAST; (ii) The sub- TeV flux depends on the surrounding 
medium density and on the fraction of thermal energy carried by the magnetic field, 6b- It ranges 
from ~ 10 -7 ergcm -2 sec -1 in the case of typical ISM density and e B < 10~ 4 to KT 10 erg cm 2 sec 1 
in the case of a source surrounded by a wind and €b ~ 10~ 05 ; (iii) The sub- TeV flux is detectable by 
high energy 7-ray experiments such as HESS, MAGIC, MILAGRO, and VERITAS; (iv) Combined 
~ 1 keV, ~ 1 GeV and sub- TeV observations will allow to determine both es and the ambient medium 
density; (v) The spectra depend only weakly on the spectral index of the energy distribution of the 
accelerated electrons. Pion production energy loss of high energy protons may contribute significantly 
in the wind case to the luminosity of high energy photons. However, it is difficult to distinguish in this 
case between the electron and proton contributions since the spectral shape is determined primarily 
by the energy dependence of the pair production optical depth. 

Subject headings: gamma rays: bursts — gamma rays: theory — radiation mechanisms: nonthermal 



1. INTRODUCTION 

In the standard fireball model of gamma-ray-bursts 
(GRB's), the observable effects are due to the dissi- 
pation of kinetic energy in a highly relativistic fire- 
ball (for reviews , see, e.g.. iPiran 11200(1 IMeszaros Il2002t 
IWaxmanl 120031) . The observed radiation is well ex- 
plained as synchrotron and inverse-Compton emission 
from shock accelerated electrons. Electrons accelerated 
in internal shock waves within the expanding fireball pro- 
duce the prompt 7-ray emission, while electrons accel- 
erated in the external shock wave driven by the fire- 
ball into the surrounding medium produce the after- 
glow emission, from the X to t he radio bands, that 
can last months after the bur st llPaczvriski fc Rhoads I 
19931 IMeszaros fc Rees"|ll997t fVietri I Il997al IWaxman I 
1997al IWaxmanl Il997ti ISari. Piran fc Naravan I 119981 
Gruzinov fc W axman ) |1999|) . 

At the early stages of the afterglow, as the en- 
ergy density of the surrounding matter swept by 
the external shock wave equals the energy density 
of the fireball plasma, a reverse shock crosses the 
fireball plasma l|Meszaros. Rees fc Panathanassioulll994t 
ISari fc Piran 1 119951) . At this transition phase, lasting 
over a duration comparable to that of the burst it- 
self, both forward and reverse shocks accelerate par- 
ticles, which subsequently emit the early afterglow. 
Early optical emission from several G RBs has been 
identi fi ed as reverse shock emission llAker lqf et al~\ 
(119991): ISari fc Piran I il999h: IMeszaros fc Rees I 1119991): 
iKobayashil II2000D: Ikobavashi fc Sari I 1200(1 ): iWei I 
(2003|): lZhang. Kobavashi fc Meszaros 1 1)2003|) ~ 
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In the past few years, evidence accumulated for the ex- 
istence of a high energy emission component (~ 1 GeV — 
1 TeV) associated with gamma-ray bursts (GRBs), which 
are typically observed at much lower energies, few 
keV to few MeV. Such a component was first ob- 
served by EGRET in the ~ 100 MeV - 18 GeV range 
(|Schneid et ol. Ill992t Ulurlev et ol. 1119941) . There is evi- 
dence for even higher energy, ~ 1 TeV emission, observed 
in one burst out of 54 BATSE GRBs in the MILAGRITO 
field of view i|Atkins et ol. II2TM l200l . 

The mechanism responsible for high energy pho- 
ton emission is still controversial. High energy 
electrons inverse Compton (IC) scatter low en- 
ergy photons to high energies, 1 GeV- 1 TeV 
ijMeszaros. Rees fc Papathanassioul|19 94) extending the 
spectrum to high energy. The exact shape of the 
spectrum depends on the values of the un c ertain 
parameters of the model jC hiang fc^JDermerJ 
iWang. Dai feLul feOQl: Zh ang fc Meszaros 1 1200 1|) . An 
additional source of high energy emission is the possi- 
ble acc el eration o f bary ons, fi rst suggested bv IWaxman I 
<|1995a|) . IVietri I l(l995|) and iMilgrom fc Usovl (|1995). 
Accelerated protons can emit synchrotron radiation 
at this energy band (IVietri I ll997bT: iTotani I Il998albl : 
IBotcher fc Dermer I Il998|) . High energy baryons can 
also produce energetic pions, via photo- meson inter- 
actions with the low energy photons, creating an al- 
ternative source of hi gh en e rgy ph otons and neutrinos 
llWaxman fc Bahcall I 119971 1200(1 IBotcher fc Dermer I 
119981) . A high energy component may also be created 
by the decay of energetic pions, produced in proton- 
proton(neutron) collisions, if the expanding plasma col- 
lides with a dense (n ~ 10 11 cm " f) cloud JKatz 1 119941: 
iDe Paolis. Ingresso fc Orlando ll200(il) . 

In this paper, we analyze the high energy emission 
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component during the early afterglow (transition phase) 
on a time scale of tens to hundreds of seconds following 
the GRB, within the framework of the fireball model. 
We explore the dependence of photon spectrum on un- 
certain model parameters, in particular on the energy 
density of the magnetic field, the power law index of 
the accelerated particles and the density of the medium 
surrounding the fireball. The ambient medium density 
differs widely between different scenarios for GRB pro- 
duction. In the neutron star merger scenar io ijGoodman I 
ll986UPaczvhski 119861 lEichler et al. I1989J) . for example, 
a density typical to the ISM (n w 1cm -3 ) is generally 
expected, while i n the massive stellar collapse scenario 
iWooslev I Il993t iLevinson fc Eichlerl Il993t iPaczvriskil 
1998) a much higher density (n « 10 3 — 10 4 cm -3 ) may 
be expected in a wind surrounding the progenitor, due 
to mass loss preceding the collapse llChevalier 1120011) . In 
addition, we explore the possible contribution of high en- 
ergy protons. 

Calculation of the early high energy emission spec- 
trum is complicated due to several reasons: First, IC 
scattering is partly in the Klein-Nishina regime. Sec- 
ond, high energy photons (produced by IC scattering 
or via pion decay) initiate high energy electro-magnetic 
cascades, the evolution of which is difficult to study 
a nalytically. Studies of cas c ade processes in the past 
( Bon ometto fc Rees I 1119711) : lOuilbert. Fabian fc Reesl 
( 1983): lSvensson I (|1987f) ^ have shown that it may have a 
significant effect on the spectrum at high energies. And 
third, the optical depth to pair production may be large 
at high energy, leading to the formation of a large num- 
ber of e^- pairs which affect the resulting spectra. Due to 
these complications, numerical calculations are required 
to accurately determine the high energy spectra and their 
dependence on uncertain model parameters. 

Numerical methods where extensively used in the past 
in the study of active galactic nucleus (AGN) plasma 
(IGuilbert. Fabian fc Rees I 1983t iZdziarski fc Lightman I 
1 19851 iFabian et al. 11198a iLightman fc Zdziarski Ill987t 
ICOTroiTIT992Tr Using these methods, new results were 
obtained, such as the complex pattern of the spectral 
indices in the X-ray (2 — lOkeV) range, that were not 
obtained by previous analytic calculations. In the con- 
text of the high energy emission from GRB's, numer- 
ical calculations are challenging. The large difference 
between the characteristic time scales for the evolution 
of high energy and low energy particles, and the fact 
that the distribution of low energy particles evolves on 
a time scale comparable to the dynamical time and 
does not reach a steady-state, implies that the numer- 
ical code must follow processes characterized by widely 
differing time scales. Due to these difficulties, the accu- 
racy of existing numerical calculations of emission from 
GRB's l)Panaitescu fc Meszaros lll998j) is limited above 
~ 1 GeV. Here we present numerical results using a code 
that allows to overcome the difficulties mentioned above. 
Our time dependent model includes an accurate descrip- 
tion of synchrotron emission, direct and inverse Compton 
scattering, pair production and the evolution of high en- 
ergy electromagnetic cascades. 

Several new high energy photon detectors will be- 
come fully operational in the near future. These in- 



clude the GLAST 3 satellite, which will greatly improve 
the sensitivity at ~ 1 GeV and will open a new win- 
dow of observation up to ~ 100 GeV. GLAST sensitiv- 
ity, ~ 10~ 13 erg cm -2 sec -1 at 1 GeV will allow detec- 
tion of many hundreds of bursts per year, if GRB's emit 
equal amount of energy at GeV and keV bands. In addi- 
tion, new generation of sub-TeV Cerenkov detectors such 
as MAGIC 4 , HESS 5 VERITAS 6 and CANGAROO III 7 
will open new era in GRB observations. The universe is 
tran sparent to 100 GeV photon s up to redshift of w 0.5 
(e.g. iSalamon fc Stecker Ill99 8[). Given the local rate of 
GRBs (|Guetta. Piran fc Waxman 1 12003) . the sensitivity 
of these detectors at ~ 100 GeV, ~ 10~ 10 erg cm -2 sec -1 
for a 100 second burst, will allow a detection of few to few 
tens of GRB's per year at this energy band (depending on 
the exact instrument energy threshold) if the GRB spec- 
trum in the keV - sub TeV range is flat, vF v tx v°. The 
field of view of the Cerenkov telescopes is small, which 
implies that fast slewing to the GRB position, on minute 
time scale, would be required. This may be achievable 
with the fast alerts that will be provided by SWIFT 8 . 
The sensitivity of the MILAGRO 9 shower detector at 
~ 1 TeV would allow it to detect ~ 1 event per year. 
The results presented here may therefore be useful for 
planning the observing strategy of the detectors and may 
allow to use high energy data to constrain the values of 
uncertain parameters of the model. 

This paper is organized as follows. In [|21we derive the 
plasma conditions during the transition phase, which is 
the duration over which the reverse shock exists. We 
calculate the critical synchrotron frequencies and lumi- 
nosities in the "wind" scenario f i!2.1|) and in the "uniform 
density ISM" scenario f H2.2J) . In i !2.3l we discuss proton 
energy loss and its contribution to high energy photon 
emission. Our numerical model is briefly presented in 
%5 A detailed description of the model may be found 
in iPe'er fc Waxman I ((2004) . In 21 we present our nu- 
merical results. We summarize and conclude in [0 with 
special emphasis on implications for high energy photon 
telescopes. 

2. MODEL ASSUMPTIONS, PLASMA CONDITIONS AND 
PROTON ENERGY LOSS AT THE TRANSITION PHASE 

The interaction of fireball ejecta with surrounding gas 
produces a reverse shock which propagates into and de- 
celerates the fireball ejecta. As the reverse shock crosses 
the ejecta, it erases the memory of the initial conditions, 
and the expansion the n approaches the Blandford-M cKee 
self-similar solutions (Blandford & McKcc 1976). We 
derive in this section the plasma parameters, the low en- 
ergy photon luminosity and spectrum, and the proton 
energy loss time scales expected during this transition 
phase, for fireball expanding into a wind and into a uni- 
form medium with density typical to the ISM. 

2.1. Wind scenario 

3 http://www-glast.stanford.edu 

4 http:/ /hegral. mppmu.mpg.de/MAGICWeb 

5 http://www.mpi-hd.mpg.de/hfm/HESS/HESS.html 

6 http:/ /veritas. sao.arizona.edu/ 

7 http:/ /icrhp9. icrr.u-tokyo.ac.jp/ 

8 http://www.swift.psu.edu/ 

9 http://www.lanl.gov/milagro/ 
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The mass loss preceding a collapse creates a density 
profile p(r) = Ar~ 2 , where the proportionality con- 
stant is taken as A = 5 x lO 11 ^* gem" 1 , and a value 
of A* ~ 1 i s assumed for a typical Wolf-Rayet star 
( Willis"! Il991l) . During the self similar expansion, the 
Lorentz factor of plasma behind the forward shock is 

Tbm = (9£ , /167re(r)) 1 ^ 2 r~ 3 / 2 , where e(r) is the energy 
density of the surrounding medium, e(r) « p(r)c 2 , and 
E is the (isotropically equivalent) fireball energy. The 
characteristic time at which radiation emitted by shocked 
plasma at ra dius r is observed by a distant observer is 
t w r/4Tl M c l|Waxman 111997c]) . 

The reverse shock is relativistic if the energy density 
of the surrounding matter swept by the forward shock 
equals the energy density of the propagating plasma. 
Since both energy densities decrease with the second 
power of the radius, a relativistic shock is formed if 
E/(Anr 2 T 2 cT) < 4T?e(r), i.e. if 

E < 167rr 2 rf cTe(r) = 5.5 x 10 56 AT 2 r 4 2 5 erg, (1) 

where T; = 10 2 - 5 I\ 2 .5 is the original ejecta Lorentz fac- 
tor, and T — 100 T 2 s is the burst duration. At higher 
explosion energy, the reverse shock is not relativistic. 
Both prompt and afterglow observations suggest that 
E < 10 54 erg, implying that a relativistic reverse shock 
should in general be formed. 

The duration T (measured in the observer frame) of the 
transition phase, during which the reverse shock exists, is 
comparable to the longer of the tw o time scales set by the 
initial conditions ( Waxman 2003J): The (observed) GRB 
duration Tqrb and the (observed) time Tr at which the 
self-similar Lorentz factor Tbm equals the original ejecta 
Lorentz factor L;, Fbm(Tt) = This implies that the 
transition radius to self similar expansion, r s = 4r BM cT, 
is 



max (2.2 x 10 16 E^T^A^ 2 : 



iAxlO^E^TT^A- 1 ) 



cm, 



(2) 



where E = 10 53 -Es3 erg. During the transition, shocked 
plasma expands with Lorentz factor close to that given 
by the self-similar solution, L ~ Fbm(?" = f s ), i.e., 



/ 9E \ 1/4 
\64irATc 3 J 



1/4^-1/4^-1/4 



(3) 



The Lorentz factor of the reverse shock in the frame of 
the unshocked plasma is T r — 1 w F^/F. 

We denote by e e and es the fractions of the thermal 
energy density that are carried, respectively, by elec- 
trons and magnetic fields. The characteristic Lorentz 
factor of electrons accelerated at the forward shock is 
Icharj = £ e rBMW p /m e . Electrons accelerated at the 
reverse shock are characterized by j c har,r = £e(r r — 
l)m p /m e ~ e e (Li/LBM)?7ip/m e . The minimum Lorentz 
factor of a power law accelerated electrons is 



"ichar X 



log (^) 1 p = 2 



(4) 



where p is the power law index of the accelerated elec- 
trons energy distribution, n e (s e ) oc s~ p . 



Synchrotron emission peaks at 
(3/2)mjl in qB/m e c. Using 



-ob. 



AttA 3 c 3 
9ET 3 



7.0 x 10 3 E 



5:5 



l/2j,-3/2 A 3/2 



A* erg cm 



-3 



( 5 ) 

for the energy density behind the forward shock at the 
transition radius, the synchrotron emission peak at the 
reverse and forward shocks is given by 



e*^ r = 0.5(1 + z)- 



7 -l/2 T -l/2 . p 2 2 1/2 v 
^53 J 2 A * 1 i.2.5 £ e.-l e B.-2 e V > 

e°* mJ = 20 (1 + z)- 1 El.i 2 T- 3/2 <-i^ 2 -2 eV, 

(6) 

where p = 2, log(7 max /7 min ) ~ 10 assumed, e e = 
0.1 e e ,-i, and es = 0.01 en, -2- These energies are above 
E. c , the characteristic synchrotron energy of electrons 
for which the synchrotron cooling time, 9m 3 c 5 / Aq i B 2 'y , 
is comparable to the ejecta (rest frame) expansion time, 
~ r/Lc, 



e° b c = 2AxlQ- 2 (l + z) 



1 El! 2 Tl' 2 A- 2 e B %^. (7) 



s° m r is comparable to the synchrotron self absorption 
energy, below which the optical depth becomes larger 
than 1, 



.ob. 



0.4(1 



T^Al'^'t^eV (8) 



(for p = 2). Equating the particle acceleration time, 
t acc ~ Ej(cqB) and the synchrotron cooling time, 
gives the maximum Lorentz factor of the accelerated 
electrons (in the plasma rest frame), 7e imax = 1-8 x 

10 7 Ell S T 2 /8 A^ /?, e B l/ _\. Synchrotron emission from 
these electrons peaks at 

e * t = 1.0xl0 u, (l + *)- 1 EU%-^A-^eV. (9) 
The reverse shock specific luminosity L e — dL/de° b ' at 

ob. ■ 



2.8 x 10 60 E^T^Al^r^eZLe-^U- 1 



-i/ 2 , 

i,2.5 c e,-l c B,-2 ' 



(10) 

For a power law energy distribution with p = 2 of the 
accelerated electrons, the specific luminosity is L £ oc e^ 1 
above e 1 ^ m .r- 

2.2. Uniform density ISM 

During self similar expansion into uniform density 
medium, the Lorentz factor of plasma is given by Tbm = 

(17E/I6imm p c?) 1/2 r~ 3 / 2 l|Blandford fc McKee I Il97^1 . 
The transition to self similar expansion occurs at a ra- 
dius 



r s = max I 6.3 x 10 b E 5 $ n 



16 pl/3„-l/3 r -2/3 



1.3xlO"^V 1/4 r a 1/4 Jcm, 



(11) 



where a typical value n = Iuq cm~ 3 of the ambient 
number density assumed. Calculations of the char- 
acteristic Lorentz factors, synchrotron frequencies and 
characteristic luminosities closely follow the steps of 
the wind scenario calculations, and partly appear in 
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IWaxman fc Bahcall I 1)2000)) . The peak energy of syn- 
chrotron emission from electrons accelerated at the re- 
verse shock, 

1/2 



£ °,m,r — 2 X 10 



(1+*)- 



r7 1/2 r 2 e 2 

n i i,2.5 e e,- 



l c B,-2 



eV, 



( 12 ) 

is comparable to the synchrotron self absorption energy, 



ef sa = 5 x 10- 2 



and is below 



(1 + *)- 



E 



1/3^1/3^-2/3^1/3 ^1/3 



53 "0 



C e,-1 C B, 



= 10 (1 + z)' 



-1/2 ^-1/2-3/2 y 
^53 ll ± 1 fc S,-2 ev - 



2 eV, 
(13) 

(14) 



Synchrotron emission from electrons accelerated at the 
forward shock peaks at higher energy, 



r ob. 



25(1 



EU 2 T 2 - e 



3/2 2 1/2 



-1 C B, 



eV. (15) 



The reverse shock specific luminosity, L e = dL/de° b 

ob. 



at 



s; c is 



L c = 1.5 x lO 59 ^ 2 T 2 1/2 



3/2 

"■0 e e,-l< -B -2 f 



(16) 



For a power law energy distribution with p = 2 of the 
accelerated electrons, L e cx e" 1 at higher energies, £ > 



2.3. Proton energy loss and its contribution to high 
energy photon emission 

Protons accelerated to high energy in the reverse shock 
contribute to the emission of high energy photons by syn- 
chrotron emission and photo- production of pions, which 
decay to produce high energy photons and electrons (and 
neutrinos). The maximum energy to which protons are 
accelerated is determined by equating the proton accel- 
eration time (in the plasma frame), t acc ~ 7 p m p c 2 / cqB , 
to the minimum of the dynamical time, the synchrotron 
energy loss time, and the pion production energy loss 
time. For expansion into a wind, the acceleration time 
is equal to the dynamical time, Td yn = 4rT, for protons 
with (plasma frame) Lorentz factor 



1p,dyn 



10 /. 1/2 1/2 



2 x 10 w A L Je 



B,-l> 



(17) 



and to the synchrotron loss time, t syn — 9m p c 5 /4q 4 B 2 j p , 
for protons with (plasma frame) Lorentz factor 



2 x 10 J 



i pl/8 T 3/8 4-3/8 

£/ 53 J 2 A* e 



1/4 
B,-l- 



(18) 



The time for e nergy loss via pion product ion may be 
approximated as (Waxman & Bahcall 1997) 



t: 



1 

2? 



f'OC 

C 

p J eo 



dea v (e)£ i (e)e 



dxx~ 2 n{x). (19) 



e/27 P 



0V (e) is the cross section for pion production for a pho- 
ton with energy e in the proton rest frame, £(e) is 
the average fraction of energy lost to the pion, and 
£o = 0.15 GeV is the threshold energy. The specific pho- 
ton density, n(x), is related to the observed luminosity 
by n(x) = L e (Tx)/4:irr 2 cTx. Photo-meson production 
is dominated by interaction with photons in the energy 
range e° b - > Teo/2j Ptdyn 



£~"mr> where L e oc e 1 l , thus 

_o6 



2 Pm^ip ^^,m,r 0~peak^peakAc 



3irr 2 T 2 e pea k 



£peak 



(20) 



A contribution from the A-resonance comparable to that 
of photons of higher energy was assumed in evaluating 
the first integral. Using a pea k — 5 x 10~ 28 cm 2 and 
^ pe ak — 0.2 at the resonance e = e vea k = 0.3 GeV, and 
Ae ~ 0.2 GeV for the peak width, the acceleration time 
equals the pion production loss time for protons with 
(plasma frame) Lorentz factor 



7p)7r « 4.6 x 10 s Eli'T^A-^Jh- 1 ^ 



Comparing eqs. (|17|l . (|18(l and 1)21)1 we find that proton 
acceleration is limited by energy losses, and that pion 
production losses dominates over synchrotron losses for 

e e > eu/10. (22) 

Note, that this is valid independent of the proton en- 
ergy, since t v oc t syn oc 7" 1 . X-ray afterglow observa- 
tions suggest that e P is close to equipartition, e s ~ 1/3 
llFreedman fc Waxman I [20011: iBerEfer. Kulkarni fc Fraill 
2003). Thus, we expect proton energy losses to be domi- 
nated by pion production. Comparing eqs. l)TT)l and lt2"T)) 
we find that for e e ~ 0.1 the highest energy protons lose 
all their energy to pion production. Assuming a power 
law distribution of proton energies, dnp/d^p oc 7~ 2 , 
roughly 0.1 of the energy carried by protons will be con- 
verted in this case to pions, and roughly half of this en- 
ergy would be converted to high energy photons through 
the decay to photons and electrons (positrons). Since 
shock accelerated electrons lose all their energy to radia- 
tion, this implies that the contribution of protons to the 
luminosity is similar to that of the electrons. This is valid 
also for values of e e well below equipartition, e e <C 0.1, 
since the contribution to the luminosity of both electrons 
and protons is oc e^ 1 in this case. 

In the case that eq. 1)22)1 is not satisfied, and the proton 
energy loss is dominated by synchrotron losses, compar- 
ing eqs. 1)17)1 and 1)18)1 implies that the highest energy 
protons lose all their energy to synchrotron losses if the 
magnetic field is close to equipartition. For lower values 

3/2 

of tB , the highest energy protons lose a fraction ~ e B of 
their energy by synchrotron emission. In this case, pro- 
ton synchrotron emission would dominate the electron 

3/2 

luminosity if e B > 10e e . 

For explosion into a uniform density ISM, similar argu- 
ments imply that pion production losses dominate over 
synchrotron losses when eq . (1221 is sat isfied . Follow- 
ing the analysis of IWaxman fc Bahcall I l)2000D . protons 
which are expected to be accelerated up to j P ,dyn ~ 

x , lose a fraction / ff (7 P ) ~ 

n ' ± 2 ,8 e e ,-i of their energy to pion 
production. Therefore, for a power law index p = 2 
of the accelerated protons, pions receive a fraction ~ 

10-3 EU S n 7 /S T} /8 ' " 1/2 



5 x IV Eli\y%'*e)il 
2xl0-S^3 /8 -o /8 T 2 - 1/8 £ 



e e ,_ie B _i of the total proton energy. 
Since in this scenario as well electrons lose almost all their 
energy to radiation, the total proton contribution to the 
photon flux is w 5 x 10~ 3 of the electron contribution, 
independent on the value of e e . If eq. 1)22)1 is not satis- 
fied, proton synchrotron losses dominate over pion pro- 
duction. A similar calculation to the former case shows 
that proton synchrotron emission dominates the electron 

3/2 a 

luminosity in this scenario only if e B > 10 e e . 

The discussion presented above demonstrates that the 
energy loss of protons is expected to be dominated by 
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pion production, and that this energy loss may produce a 
luminosity similar to that of the electrons in the scenario 
of explosion into a wind. Finally, it should be pointed 
out that if protons are accelerated to a power-law energy 
distribution with an index p significantly larger than 2, 
their contribution to the luminosity will be significantly 
reduced, since the fraction of energy carried by the high- 
est energy protons will be small. 

3. THE NUMERICAL MODEL 

The acceleration of particles in the two shock waves 
that exist during the transition phase is accompanied 
by numerous radiative processes. In the numerical cal- 
culations, w e use the time dependen t numerical model 
described in lPe'er Sz Waxman I l)2004f) . Our time depen- 
dent model follows the evolution of the particle distri- 
bution and the emergent spectra, by solving the kinetic 
equations for the electron and photon distributions, tak- 
ing into account synchrotron emission, inverse Compton 
scattering, pair production and pion photo-production 
interactions. These calculations are done for a wide 
range of particle energies, including the evolution of rapid 
electro-magnetic cascades at high energies. 

In the calculations, we focus on particles that pass 
through the forward or the reverse shock waves and 
therefore are in the downstream region of the flow, rela- 
tive to the relevant shock wave. While the upstream rel- 
ativistic flow is highly anisotropic, the shocked gas ther- 
malize, hence isotropise (in the plasma comoving frame) 
immediately after passing the shock, on a characteristic 
length scale of several skin-dep th fe.g.. lKirk et al. Ill998l 
l2000HFrederiksen et al. 112004^ ■ The emitted radiation is 
therefore isotropic in the comoving frame. Most of the 
shocked gas and most of the blast-wave energy are con- 
centrated in a shell of co-moving thickness Ar = Cc Xrf, m 
with C, < 1 (e.g., in the lBlandford fc McKeel l)1976j) self- 
similar solutions which give the spatial dependence of 
the hydrodynamical variables, 90% of the energy is con- 
centrated in a shell of co-moving thickness correspond- 
ing to C = 1/7). Since the details of the spatial depen- 
dence of the electron and magnetic field energy fractions 
are not known, we adopt in our model the commonly 
used approx i mation (e.g.. iLightman fe; Zdziarski I Il987t 
iConpj 1119921 iPilla fc Loeb 111998ft that radiation is pro- 
duced within a homogeneous shell of comoving width 
Ar = ^cTdyn with £ = 1. Since the shocks velocity is 
time independent during Td yn , the shock- heated comov- 
ing plasma volume is assumed to increase linearly with 
tim e, i.e., constant particle number density is assumed 
fsee lPe'er fc Waxman"ll2004 for further details). 

In the comoving frame, homogeneous and isotropic dis- 
tributions of both particles and photons are therefore as- 
sumed. Parallel calculations of processes in the two shock 
waves are carried out, where photons produced at each 
shock participate in IC scattering, photo-production in- 
teractions and pair production interactions occurring at 
both shocks. 

The particle distributions are discretized, the proton 
spectrum spans 11 decades of energy < 10 11 ), and 
the electron spectrum 14 decades (j e < 10 14 ). A fixed 
time step is chosen, typically 10~ 4 times the dynamical 
time. Numerical integration is carried out with this fixed 
time step. Particles and photons for which the energy 
loss time or annihilation time are shorter than the fixed 



time step, are assumed to lose all their energy in a single 
time step, producing secondaries which are treated as a 
source of lower energy particles. Photo-pion production 
is calculated at each time step by direct integration of 
the second integral in Eq. 1191 while approximating the 
first integral by the contribution from the A-resonancc 
(see the discussion following Eq. I2U|) . Half of the energy 
lost by protons goes into ir° that decays into 2 photons, 
each carrying 10% of the initial proton energy. 

The following approximations are made: (i) Plasma 
parameters are assumed to be time independent during 
the transition phase; (ii) The fraction of thermal energy 
that is carried by electrons (magnetic field), e e (es), is 
the same at the forward and reverse shock waves; (iii) 
The power law index of the energy distribution of accel- 
erated electrons and protons is the same at both shocks. 

In the present calculations, we do not consider proton 
synchrotron emission, which is dominated by pion pro- 
duction for the considered parameter range. Synchrotron 
self absorption is also not considered here, being irrele- 
vant for processes occurring at high energies. Photons 
below the self absorption frequency, e° b a < 1 eV (in both 
scenarios; see eqs. 03 I13|l can produce pairs only with 
photons energetic than T 2 {m e c 2 ) 2 / 'e° b a ~ 10 16 eV, which 
are absent (see the numerical results below). More- 
over, these photons are below the threshold energy for 
pion production, e^g > 0.2rGeV/7 P = 10 eV, where 
j Pt7t — 5 x 10 9 (see eq. |?TJ, and T = 10 2 - 5 where taken. 

4. NUMERICAL RESULTS 
4.1. High vs. low density 

Figure^shows spectra obtained from a GRB transition 
phase in the two scenarios discussed: Explosion in a uni- 
form low density medium typical to the ISM, n = 1 cm -3 
, and explosion in a A* — 1 wind, where the density of 
ambient medium at the transition radius is much higher, 
n ~ 10 3 cm -3 . 

The most evident feature of the two spectra is their 
flatness, vF v oc v a with a ~ 0, extending over seven 
energy decades between e° b m j ~ 1.4 keV and £° b max — 
60 GeV for expansion into ISM, and over ten energy 
decades from e° b mr « 1 eV to £° b max ~ 10 GeV for explo- 
sion into a wind. This is a result of dominant synchrotron 
emission term at these energies. 

The main difference between the two scenarios appears 
at high energies, 100 GeV — 1 TeV, due to the different 
optical depth to pair production, r 77 (e), as shown in Fig- 
ure |3 The large optical depth to pair production in the 
wind case softens the high energy IC spectrum, result- 
ing in spectral index a = —1. Consequently, while in 
the ISM case the 1 TeV flux is comparable to the flux at 
1 GeV, in the high density wind case the 1 TeV flux is 
2.5 orders of magnitude lower. A second difference be- 
tween the two scenarios appears at the optical-UV band. 
In the low density scenario, at low energies, e < s° m f, 
emission from both shock waves produces a moderate in- 
crease in the spectral slope, a ~ 0.3 above e° b c « 20 eV. 
As demonstrated in Figure this slope results from a 
flat spectrum produced at the reverse shock, a — above 
£ ° b m,r ~ 10 -2 e V, and forward shock emission character- 
ized by a — 0.5 below e° b m ^. In a high density medium, 
a nearly flat spectrum is obtained above e ob m r ~ 1 eV. 
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In an explosion into constant density ISM scenario, the 
total energy carried by the IC scattered photons is w 10% 
of the total energy in synchrotron emitted photons, well 
below the simple estimate (ce/es) 1 ^ 2 , which is based on 
the assumption that all the synchrotron emitted pho- 
tons serve as seed photons for Compton scattering. This 
discrepancy is due to Klein-Nishina suppression of the 
cross-section at high energies. Consequently, IC scatter- 
ing becomes dominant only above ~ 50 GeV, where a flat 
spectrum is created for a power law index p — 2 of the 
accelerated electrons energy distribution. In the wind 
case, IC scattering is more prominent due to the higher 
photon density. 

Lightcurve of the spectrum of an explosion into con- 
stant density ISM scenario is presented in Figure 0] The 
constant injection rate of energetic electrons results in 
a linear increase of the synchrotron emission component 
in time. Since the energy of photons emitted by elec- 
trons for which the synchrotron cooling time equals the 
time t, £ 7 . c (i) oc t~ 2 , the energy density in the low en- 
ergy photons increases with time. Only photons with 
e° < 1 keV are in the Thompson limit for Compton 
scattering, therefore the inverse Compton component in- 
creases in time, at early times. At later times, t < Td yn , 
where Td yn — ATbmT is the dynamical time of the prob- 
lem, the low energy photon density becomes high enough 
that pair production phenomenon limits the high energy 
emission component. 

The contribution of protons to the photon flux is neg- 
ligible in the ISM case, however this contribution is 85% 
of the electrons' contribution to the flux in the wind 
scenario (see H2.3JI . Most of this contribution results 
from the decay of energetic pions, initiating high en- 
ergy electro-magnetic cascades. The high energy cas- 
cade results in flat (a — 0) spectrum at the observed 
energy range, which, for power law index of the accel- 
erated particles p — 2 is similar to the spectrum from 
inverse Compton scattering. An additional small contri- 
bution from the decay of low energy pions results in a 
spectral slope a = 1 , contributing to the small deviation 
from flat spectrum shown below 1 GeV in this scenario. 
Above the energy at which the optical depth for pair pro- 
duction r 77 > 1, ~ 1 GeV, the spectrum is determined 
by the spectrum of low energy photons, thus the proton 
contribution is hard to discriminate at these energies as 
well. 

4.2. Fraction of energy carried by electrons and 
magnetic field 

In an explosion into a low density ISM, 7 TOjr < 7 C <C 
7max- Thus, for p ~ 2 electrons lose nearly 100% of their 
energy to radiation, resulting in a nearly linear depen- 
dence of the flux on e e . The dependence of the spec- 
trum on the fraction of thermal energy carried by the 
magnetic field, es, is illustrated in Figure The dom- 
ination of synchrotron emission implies that the values 
of e 7jC and e T)m determine the spectral behavior at low 
energies, while a flat spectrum is expected for p = 2 at 
higher energies, up to £ 7 . max ~ 50 GeV. The ratio of 
1 GeV to 1 keV flux is informative about the value of cb , 
with large ratio (~ 30) implying small e#(~ 10~ 4 ). At 
higher energies, e > 100 GeV, a flux comparable to the 
1 GeV flux is produced by IC scattering for es < 0.1. 

For explosion into a wind, 7 C < j m . r *C 7max- Here too, 
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Fig. 1. — Predicted photon spectrum during the transition phase. 
Results are shown for E = 3 X 10 53 erg , T; = 10 25 , T = 10s, 
e e = 10 -1 , £b = 10 — 2 , p = 2.0. Solid: explosion into uniform 
low density medium (ISM), n = 1cm -3 . Dashed: Explosion into 
a wind with A, = lgrcm - 1 . Dotted: Explosion into a wind with 
A» = lgrcm -1 , proton contribution to the flux omitted. Lumi- 
nosity distance = 2 X 10 28 cm and z = 1 were assumed. Inter- 
galactic absorption of high energy photons through pair production 
interactions with the IR, background, which becomes appreciable 
above ~ 0.1 TeV, is not taken into account. 




Obs. energy [eV] 

Fig. 2. — Energy dependent optical depth to pair production, 
for the two scenarios considered in figure Solid: explosion into 
ISM, dashed: explosion into a wind. All physical parameters are 
the same as in figure HI 




Fig. 3. — Two emission components of the observed spectra, ex- 
plosion into ISM. Dash: forward shock emission, dash-dotted: re- 
verse shock emission, solid: combined spectra. Left: equipartition, 
e e = 10 -0 ' 5 , £b = 10 -0 ' 5 , p = 2.0. The two emission components 
are comparable. Right: e e = 10 -1 , eg = 10 -2 , p = 3.0. The wavy 
shape results from different peaks of the two synchrotron compo- 
nents, and the two IC components. All other physical parameters 
are the same as in figure HI 



for p ~ 2 electrons lose nearly 100% of their energy to ra- 
diation, resulting in a nearly linear dependence of flux on 
e e . The dependence of the spectrum on eg is illustrated 
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Fig. 4. — Lightcurve of the spectra for the scenario of explosion 
into uniform low density medium (same parameters as in figure 0. 
The lines show the spectrum after fractions of 1/100, 1/30, 1/10, 
1/3 and 1 dynamical time. 
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Fig. 6. — Dependence of the spectra on the fraction of thermal 
energy carried by the electrons and by the magnetic field, e e and 
tB, explosion into a wind characterized by A* = lgrcm - 1 . Solid: 
e e = 10~ ' 5 , e B = 10~ ' 5 , dashed: e e = 10" 1 , e B = 10~ 2 , dash- 
dotted: t e = 10 _1 , eg = 10 -4 . All other physical parameters are 
the same as in figure 111 



in Figure [fj] Here, IC emission is much more prominent. 
IC energy loss of the electrons affects the spectrum at all 
energy bands by modifying the electrons energy distri- 
bution. The spectrum becomes harder with lower values 
of €b, reaching a ~ 0.2 below lOGeV for es = 10~ 4 . In 
this scenario, the minimum Lorentz factor of electrons 
accelerated at the front shock, 7 m i n ./ = 2.0 x 10 3 (see 
Eq. is comparable to jic-s ~ 10 4 , which is defined 
as the Lorentz factor of electrons that emit synchrotron 
radiation above their Klein-Nishina limit. We show in 
Appendix [S] that in this case, if the main energy loss 
mechanism is inverse Compton scattering, a spectral in- 
dex in the range < a < 2/3 should be expected. A 
1 GeV to 1 keV flux ratio larger than 10 indicates there- 
fore e_B ~ 10 -4 , and smaller ratio indicates higher value 
of €b, similar to the case of explosion into low density 
ISM. Above ~ 10 GeV, the flux is suppressed by pair 
production. 




10 10 10 
Obs. energy [eV] 

Fig. 5. — Dependence of the spectra on the fraction of thermal 
energy carried by the electrons and by the magnetic field, e e and 
eg, explosion into a low density ISM, n = 1cm -3 . Solid: e e = 
10 -0 ' 5 , eg = lO -0 ' 5 , dashed: e e = 10 _1 , eg = 10" 2 , dash-dotted: 
e e = 10 — 1 , eg = 10 -4 . All other physical parameters are the same 
as in figure 111 



4.3. Power law index of the accelerated particles 



The dependence of the spectrum on p varies with the 
external density, as demonstrated in Figures and [SJ 
The dominant synchrotron component in explosion into 
a low density ISM results, for large p, in a complex wavy 
shape shown in Figure [7| The four peaks result from 
synchrotron and IC components at the two shock waves, 
as demonstrated in Figure (right). The dependence on 
p is, however, not prominent. For explosion into a wind, 
the prominent IC component leads to a spectrum which 
is in nearly independent on p. 




10' 10" 10° 
Obs. energy [eV] 

Fig. 7. — Dependence of the spectra on the power law index p 
of the accelerated electrons, explosion into low density ISM, n = 
1 cm -3 . Solid: p = 2.0, dashed: p = 2.5, dash-dotted^p = 3.0. All 
other physical parameters are the same as in figure HI The wavy 
shape of the spectrum results from synchrotron and IC components 
of the two shock waves, as seen in figure O 



5. DISCUSSION 

We have presented numerical results of calculations of 
the early afterglow GRB spectra (figures Q 1518(1 within 
the fireball model framework. Our time dependent nu- 
merical code describes synchrotron emission, inverse- 
Compton scattering, e pair production, photo-pion pro- 
duction and the resulting high energy cascades. We 
have shown that the dependence of the spectra on the 
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Fig. 8. — Dependence of the spectra on the power law index 
p of the accelerated electrons, explosion into a wind with A» = 
lgrcm -1 . Solid: p = 2.0, dashed: p = 2.5, dash-dotted: p = 3.0. 
All other physical parameters are the same as in figure HI 



ambient density and the magnetic field equipartition 
fraction is pronounced mainly at high photon energies, 
1 GeV — 1 TeV (see Figures El : Comparable flux at 
1 GeV and 1 TeV implies low density and es — 1CP 4 , 
while a large ratio of ~ 10 3 results from dense medium 
and 6b near equipartition. A second indication for 
€b — 1CT 4 is a large (~ 30) ratio of the flux at 1 GeV to 
the flux at 1 keV. In all cases, the spectra depend only 
weakly on the power-law index p of the accelerated elec- 
trons. The flux expected in the 1 — 50 GeV energy band, 
vF v ~ 10 -8 — 10~ 7 ergcm~ 2 sec^ 1 for z = 1 bursts, is 
well within the detection capability of GLAST. 

The recent detection by MILAGR ITO of 1 TeV photon 
flux associated with GRB970417a ijAtkins et al. 112003ft 
can be explained as due to the onset of fireball decel- 
eration by a medium with density typical to the ISM, 
provided that the magnetic field is well below equiparti- 
tion, es < 10~ 4 . The flux at 1 TeV is expected in this 
case to be an order of magnitude higher than the 100 keV 
(BATSE) flux, as inferred from the MILAGRITO de- 
tection. We note, that an alternative explanation was 
suggested for the 1 TeV flu x, synchrotron e mission from 
shock accelerated protons ijTotani l|1998albf l. This expla- 
nation requires a very low fraction of the energy to be 



carried by electrons, e e ~ 10 -3 , and an isotropic equiva- 
lent explosion energy E ~ 10 56 erg. We find this expla- 
nation less attractive, since X-ray afterglow observations 
imply e e values near equipartition and isotropic equiva- 
lent energies in the range of ~ 10 515 erg to 1Q 53 - 5 erg 
llFreedman fc Waxman 1120011: iBerger. Kulkarni fc Fraill 
2003). 

The energy loss of high energy protons is dominated by 
photo-production of pions, rather than by synchrotron 
emission, for e e > es/10. The contribution of pion decay 
to the high energy photon luminosity may be comparable 
to that of inverse-Compton emission of shock accelerated 
electrons in the case of fireball expansion into a high den- 
sity wind, and power law index p ~ 2 of the accelerated 
particles. However, it is difficult to distinguish in this 
case between the electron and proton contributions since 
the spectral shape is determined primarily by the energy 
dependence of the pair production optical depth. In the 
two scenarios considered, explosion into constant density 
ISM and into a pre-ejected wind, if e e < 10~ D — 10~ 4 , pro- 
ton synchrotron emission becomes the dominant emission 
mechanism. 

Although the MILAGRITO detection is singular, the 
new generation of high energy sub- TeV detectors, such 
as MILAGRO, VERITAS, HEGRA, HESS, and MAGIC 
that will become fully operational in the next few years, 
will be sensitive enough to observe the high energy early 
afterglow emission at 100 GeV from ~ 10 bursts per 
year. The detection rate of higher energy, 1 TeV com- 
ponent is much lower, due to the large pair production 
optic al depth for TeV photons originating at z > 0.3 
Ce.g.. lSalamon fc Stecker ill 998(1 . However, the predicted 
1 TeV flux implies that w 1 burst per year, that arrives 
from low redshift z < 0.3, could be detected in the 1 TeV 
band as well. Such detections will allow to constrain two 
of the uncertain parameters of the fireball model: the 
ambient gas density, and the magnetic field equipartition 
fraction. 
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APPENDIX 

SPECTRAL DISTRIBUTION RESULTING FROM SYNCHROTRON RADIATION AND IC 

SCATTERING 

If inverse Compton scattering is the dominant energy loss mechanism of electrons that were accelerated with a power 
law energy distribution dn e /dy oc "f~ p , the resulting power law index of energetic electrons above max(7 m , j c ) deviates 
from p + 1 . 

The resulting power law can be calculated analytically if the minimum Lorentz factor of the accelerated electrons 
7 m , is comparable to 7/c-Sj which is defined to be the Lorentz factor of electrons that emit synchrotron radiation at 
the Klein-Nishina limit, h(3qBj] Cs ) / (2m e c) = 3/i(m e c 2 )/jic s - At steady state, a broken power law is expected with 
index (pi + 1) above Jic—S an d (P2 + 1) at lower energies. Synchrotron emission has a broken power law spectrum: 

edn 7 /de oc £ - pi/2 -» U y oc e 1 '^ 1 ' 2 (e > e W -s), 

edn 7 /de oc e~ P2/2 -> U 7 oc e 1 ^ 2 ' 2 (s<e IC -s), (Al) 
where eic-s = ^(SqB'yjQ )/(2m e c). Energetic electrons therefore lose energy at a rate dj/dt oc j 2 Uj oc <-y 1+ P 2 / 2 . As 
these electrons are injected with power law index p, their steady clist.ribut.iou is dn e (^ > yic-s)/d') 

Pi=P+7T-l- (A2) 
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At lower energies, electrons lose energy by scattering photons above Eic-s, thus dn e (^ < 7/c-s)/^7 oa 7 ( 1+ pi/ 2 ) ; or 

Pi 



P2 



(A3) 



Combined together, one obtains 



Pi = gfa - 1), 

P2 = 



(A4) 



As an example, assuming power law index p = 2 of the accelerated electrons, the photons spectral index is 1/3 at low 
energies, or vF v cx j/ 2 / 3 . 

Note that the above analysis is valid provided both p\,p2 < 2, which is translated via eq. I|A4J) into the demand 
p<2.5. 
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